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We study cosmological perturbations in Hofava-Lifshitz Gravity, a recently proposed potentially ultraviolet- 
complete quantum theory of gravity. We consider scalar metric fluctuations about a homogeneous and isotropic 
space-time. Starting from the most general metric, we work out the complete second order action for the per- 
turbations. We then make use of the residual gauge invariance and of the constraint equations to reduce the 
number of dynamical degrees of freedom. At first glance, it appears that there is an extra scalar metric degree 
of freedom. However, introducing the Sasaki-Mukhanov variable, the combination of spatial metric fluctua- 
tion and matter inhomogeneity for which the action in General Relativity has canonical form, we find that this 
variable has the standard time derivative term in the second order action, and that the extra degree of freedom 
is non-dynamical. The limit A ^ 1 is well-behaved, unlike what is obtained when performing incomplete 
analyses of cosmological fluctuations. Thus, there is no strong coupling problem for Hofava-Lifshitz gravity 
when considering cosmological solutions. We also compute the spectrum of cosmological perturbations. If the 
potential in the action is taken to be of "detailed balance" form, we find a cancelation of the highest deriva- 
tive terms in the action for the curvature fluctuations. As a consequence, the initial spectrum of perturbations 
will not be scale-invariant in a general spacetime background, in contrast to what happens when considering 
Hoi^ava-Lifshitz matter leaving the gravitational sector unperturbed. However, if we break the detailed balance 
condition, then the initial spectrum of curvature fluctuations is indeed scale-invariant on ultraviolet scales. As 
an application, we consider fluctuations in an inflationary background and draw connections with the "trans- 
Planckian problem" for cosmological perturbations. In the special case in which the potential term in the action 
is of detailed balance form and in which A = 1, the equation of motion for cosmological perturbations in the far 
UV takes the same form as in GR. However, in general the equation of motion is characterized by a modified 
dispersion relation. 
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I. INTRODUCTION 

Recently, Hofava proposed ll], 0] (see also IH) a model for quantum gravity wiiich is power-counting renormalizable and 
hence potentially ultra-violet (UV) complete. The model is based on a scaling symmetry which treats space and time differently. 
Hence, the model explicitly breaks Lorentz (and hence also general coordinate) invariance. The action chosen by Hofava is 
power-counting renormalizable with respect to the scaling symmetry ' . For a particular value of one of the coefficients in the 
Lagrangian (A = 1), the infrared (IR) limit of the action reduces to that of Einstein gravity ^. 

Since it was proposed Hofava-Lifshitz (HL) gravity, as this model is now called, has attracted a lot of attention (for a complete 
list of references, the reader is referred in the recent paper We will only mention some papers relevant to our study. Initially, 
gravitational wave solutions in HL cosmology were studied 131. The first papers on the early universe cosmology of HL gravity 
are lH, where it was realized that the analogs of the Friedmann equations in HL gravity include a term which scales as dark 
radiation and contributes negatively to the energy density. Thus, it is possible in principle to obtain a nonsingular cosmological 
evolution with the Big Bang of Standard and Inflationary Cosmology replaced by a bounce. In this context, it becomes possible 
ifioll to provide a realization of the "matter bounce" alternative ifllll to cosmological inflation for explaining the origin of an 
almost scale-invariant spectrum of cosmological perturbations. As realized originally in 10] the different ultraviolet behavior of 
the theory might provide an alternative to cosmological inflation for solving the problems of Standard Cosmology such as the 
horizon and flatness problems i^. The specific UV scaling of HL gravity could change the usual arguments for the origin of the 
scale-invariance of cosmological perturbations in inflationary cosmology, as pointed out in In fact, the dominant UV terms 
in the action may lead to the possibility of obtaining a scale-invariant spectrum of cosmological perturbations in the expanding 
phase of HL cosmology without inflation lfl2ll (see also lfl3ll ). To substantiate these conclusions, however, a careful analysis of 
the theory of cosmological perturbations in HL gravity is required, and this is the topic of the current paper. 

A possibly more important reason to study cosmological perturbations in HL gravity are basic consistency issues of HL gravity 
itself ^ . For a particular value of one of the parameters in the HL action, namely A = 1, the theory has its IR fixed point the action 
of GR. HL gravity has the same dynamical degrees of freedom as General Relativity (GR), but it does not have the complete 
diffeomorphism invariance of GR. Spatial diffeomorphisms are still a symmetry, but space-dependent time reparameterizations 
are no longer allowed. Thus, one loses one out of the four gauge modes of GR, and hence one extra physical mode is expected 
survive. This fact was already pointed out in [Ij and has been further discussed in |[^ [3, [IH [l^ . This extra physical mode, if 
dynamical, would lead to severe problems for HL gravity, since no effects of extra gravitational deg rees of freedom have been 
observed, and since there are stringent limits on the presence of such degrees of freedom. In lfl4lll6ll . cosmological fluctuations 
in the absence of matter were considered. A new physical scalar gravitational mode was found. In the limit A = 1, this mode 
was claimed to be non-dynamical 13, although the constraint equation in both lfl4ll and lfl6ll showed a singularity in this limit. 
Perturbations in the presence of matter were very recently considered in UtIi . where it was claimed that the extra gravitational 
degree of freedom is physical and becomes strongly coupled in the limit A = 1 in which GR is to be recovered. 

These differing claims form our second motivation to perform a careful study of cosmological perturbations in the presence of 
matter in HL gravity. It is crucial to perturb about a dynamical Friedmann universe as opposed to perturbing about Minkowski 
space-time Our work builds on the paper ifTsIl in which the groundwork for the present study was provided. We start with the 
general metric including scalar metric fluctuations. We use the spatial diffeomorphism invariance to choose a gauge in which 
the spatial metric is diagonal (we focus on the scalar metric fluctuations ^). Making use of the perturbed constraint equations, 
we determine the action for the remaining degrees of freedom of the cosmological perturbations. At this stage there are indeed 
two app arent physical degrees of freedom present. After expressing the action in terms of the usual Sasaki-Mukhanov variable 
ll20l I2III . the variable in terms of which the action for cosmological perturbations has canonical kinetic term, we find that the 
extra degree of freedom for scalar metric fluctuations is non-dynamical. In particular, there is no strong coupling problem for the 
fluctuations: consistently including cosmological expansion regulates the divergence found in fl%. In the case of cosmological 



' See also Q| for a recent study of the renormalizability issue. 
- See also Ql for a study of the IR limit. 

' One of the authors (RB) thanks G. Dvali and E. Witten for emphasizing this issue, and K. Zarembo for interesting discussions on this point. 

"* Doing the latter is inconsistent with the background constraint equations since the presence of any matter will lead to a non-vanishing average energy density 

and hence to cosmological expansion. 
^ See fl^ for a review of the theory of cosmological perturbations and fl9ll for a shorter overview. 
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perturbations, the absence of new dynamical degrees of freedom holds for any value of A, whereas for gravitational waves our 
analysis only covers the case A = 1, the most interesting case. 

The action for cosmological perturbations derived in this paper allows us to study the spectrum of curvature fluctuations in 
HL gravity. If the potential term in the Hofava action is taken to be of "detailed balance" form, we find a cancellation of the 
leading UV terms in the action. Thus, unlike what happens in the case of HL matter on a fixed background, the initial spectrum of 
fluctuations is not scale-independent. Scale invariance of the fluctuations in the UV region is maintained if we add terms which 
break the detailed balance condition. More generally, we study ways of obtaining a scale-invariant spectrum during the course 
of cosmological evolution. 

The outline of this paper is as follows: We begin with a brief review of HL gravity. In Section 3 we discuss cosmological 
perturbations and show that no strong coupling problem arises. In Section 4, we compute the power spectrum of cosmological 
perturbations and discuss applications to inflationary cosmology. The final section contains a discussion and conclusions. 



II. SETUP 



A. Brief review of Hof ava-Lifshitz tlieory 



The dynamical degrees of freedom in HL gravity are the usual metric degrees of freedom which appear in the ADM approach 
to canonical gravity, namely the spatial metric gij, the lapse function N and the shift vector Ni. In terms of these fields, the full 
space-time metric is 



ds^ = -N'^dt^ + g^j{dx' + N'dt) [dx^ + N'dt) , 



(2.1) 



where the indices of iV^ are raised and lowered using the spatial metric . 

Note that here we will allow N = N{t, x), as the most general form in the ADM decomposition. Because if one assumes 
N = N{t), then the Hamiltonian constraint for perturbations will be lost, and one can not recover GR in the IR limit and in the 
case A = 1. For studies of cosmological perturbations assuming N — N{t), the reader is referred to ifTil [Tsl [l6ll . 

The action of Horava-Lifshitz gravity contains a "kinetic" part and a "potential" part. 



with 



where 



qg 



Q9 — G9 I C9 



dtd^Xy/gN {K,jK'^ - XK^ 



(2.2) 



(2.3) 



is the extrinsic curvature and K = g^^ Kij. Hofava chose the potential to be of the "detailed-balance" form ' 



5^ = / dtd^x^N 



.-2 
2^2 



where is the Cotton tentsor defined by 



Akl 



8(1 -3A) 



1 -4A 



-R' 



AR - 3A' 



= — Vfc ( R] - -R5', 
V9 



(2.4) 



(2.5) 



Note that in ( 12.3b A is a dimensionless coupling of the theory and therefore runs as a function of energy. The extra terms in the 
potential contain two further constants w and ^. For our applications to cosmology, the terms involving w will not play a role. 
The action of GR is recovered in the IR limit if A = 1. 

The general structure of the action of scalar field matter in Hofava-Lifshitz gravity contains two parts; a quadratic kinetic term 
invariant under foliation-preserving diffeomorphisms and a potential term: 



5"^ = / dtd^x^N 



1 



2N^ 



{if - N'diip) + diip, g,j) 



(2.6) 



* It is argued in (^ [TtIi that solutions of GR are often not recovered if the potential term is taken to be of detailed balance form, a problem already encountered 
in (23l in the context of spherically symmetric metrics. This problem is due to a strong coupling signature which does not arise if renormahzable terms which 
break the detailed balance condition are added to the potential. Going beyond the detailed balance form of the potential also allows IR solutions with a positive 
or vanishing cosmological constant I23ll24ll . whereas maintaining the detailed balance condition yields a negative cosmological constant. 
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with the "potential terms" 

F{f,di(p,gij) = -~V{(p) + .91C1 + 511^1 + + 526 + 312^16 +536 , (2-7) 
where the £,i are invariants built out of spatial gradients of (p: 

ei = 9V5.</^, (2.8) 

6 = (A(p)2, (2.9) 

6 = (Av?)(AV). (2.10) 

Here and in what follows, we use A = d^di = i9^/a^ as a shorthand. Note that gi must be negative in order to obtain the 
standard form of the kinetic term in the IR. 



B. Background equations of motion 

The equations of motion for and A'^; are the energy constraint and momentum constraints, respectively. They take the general 
form 



i?2 + Ai? - 3AM - - N^d,ipy + F , (2.11) 



8(1 - 3A) V 4 J 2m 

= ^V, (kI - XKSj) - ^ - N^d.,ip) . 

In this work, we focus on a spatially flat background. Thus the background values for the metric are 

N=l, iV, = 0, g,^=a^5^j, ipo--Mt), (2.12) 

where a — a{t) is the traditional scale-factor In this background, one has Cij = and e^3^RaVjR\ = 0. Hence, at the 
background level, the energy constraint gives 

°-8(3A-l)~^°"T + ' ^^-^^^ 

while the momentum constraint is trivially satisfied. 

The space diagonal component of the generahzed Einstein equation takes the form 



2(3A-1) / • , ^ .X , 3a«^AV , 1 . 



where H = a/ais the Hubble parameter and Vq = V{(po). 

Combining Eq. ( 12.14b and the energy constraint ( 12.131 ) equations, we find another useful equation: 

The background equation of motion for the scalar field is as usual 

<fQ + 3H^o + V' :^0. (2.16) 



III. PERTURBATION THEORY AND THE DYNAMICAL DEGREES OF FREEDOM 



The scalar metric fluctuations about our background can be written as (see 11 1811 for an overview of the theory of cosmological 
perturbations) 

5500 = (3.1) 

Sgo^ = a^d,B, (3.2) 

= -2a\^d,j-d,djE) (3.3) 

where cj), B, ip and E are functions of space and time. Matter is perturbed, as well. The scalar field perturbation is denoted by 

Sip = Q. 
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The theory is invariant under spatial diffeomorphisms 

^ x' + f\x\t) (3.4) 

and under space-independent time reparametrizations. Compared to the situation in GR, one has lost the invariance under space- 
dependent changes in time. 

We can use the spatial diffeomorphism invariance to choose the gauge E = {], but one then no longer has the extra gauge 
freedom to choose spatially flat gauge {ip = Qm addition to E — Q) or longitudinal gauge {B — Qin addition to E — 0). We are 
left with three variables, namely 0, ijj and B. The lapse and shift functions become 

N = 1 + , 

(3.5) 

N, ^ d^B{t,x'), 

and that cj), B and ip are of the same order as the scalar field perturbation Q = Sip. As in GR, in order to get the second order 
(and even the third order) action, we only need to expand N and Ni to first order. Higher order contributions of N and Ni can 
be eliminated making use of the equations of motion. 

One should also note that in the E = gauge, the spatial metric gtj is conformal flat. In this gauge, making use of the local 
Weyl transformation, one has Cij = and e^^^ RuV jR\ = 0. So the parameter uj does not enter the cosmic perturbation theory. 

A. Constraints 

At first order, the energy constraint gives 

= 2(1 - 3A) {K^ifl + 12(1 - 3X)H^) <j) + 8(3A - 1)^HAB 

+ 2k2(3A - 1) i^ipoQ + V'Qj + kV^A— ^ + 24(3A - l)^HiP 

and the momentum constraint gives 

= (3A - l)H(l> + (A - 1)AB + (3A - 1)^' - —ifioQ . (3.7) 
From i3.6i and (13.7b we can solve for cj) and B explicitly to get 

^ = 16i/^(l-3A)^ + 2^^(A-l)(3A-l)^g {^^^'^-^ + - ^^^^^ - ''^'"^ 

+2k\-1 + A)(-l + 3A)Q<^o + 2Qk\-1 + 3A) {H{-1 + 3A)^o + (-1 + A)l/')} , 

32ff^(3A-l)+4.^(A-l)^g {-2i/A^.^A,^ - 4.^(-l + 3A)^o {hQ + ^^o) 
+Qk^ (I2ij2(l - 3X)ipo + K^^l + 4i/(l - 3A)y') } . 



(3.8) 



AB 



Note that these constraint equations allow us to solve for (p and B without any singularities, even in the case A = 1. If one 
(incorrectly) had expanded about Minkowski space-time instead of about an expanding universe, one would have obtained a 
singularity in the constraint equation (see e.g. Eq. (68) of ITtII ). Note also that and B remain perturbatively small. We learn the 
important lesson that the expansion of space which is inevitable in the presence of matter removes the potential strong coupling 
problem for cosmological perturbations. This is an important consistency check for the cosmology of HL gravity. 



B. Second-order action 



At this stage, we have two independent degrees of freedom for scalar metric fluctuations (instead of only one as would be the 
case in GR), namely ip and Q. In order to set up the linear theory of cosmological perturbations, we need to find the second order 
action for the fluctuations. 

After inserting the perturbed metric and perturbed matter into the action for gravity and matter, expanding to second order in 
the fluctuations, and make use of the constraints ( |3.8l l, we obtain 

S2[i': Q]^ j dt(fx {c^ + UQQ- giQAQ + .92 (AQf + g3{AQ)iA^Q) + 
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with gi = gi{(fo). The expressions for the various coefficients can be found in Appendix I A II 

It appears that there are two dynamical degrees of freedom, ijj and Q respectively. However, it will be shown that this is an 
illusion (see the next subsection for a detailed discussion on the issue of dynamical degrees of freedom in Hofava theory). In 
fact, there is only one dynamical degrees of freedom, the same as in GR. The easiest way to see this is to look at the combination 
of all of the kinetic terms in the above action, i.e. the Q^, ilP' and Qil> terms (see Appendix IA 11 1, and to realize that they can be 
brought into a "perfect square" form 

/ • H ^ 

\ 

This fact implies that there is indeed only one dynamical degrees of freedom in our system. This degree of freedom is precisely 
the Sasaki-Mukhanov combination of matter and metric fluctuations ll20l I21I1 . the variable in terms of which the action for 
cosmological perturbations in GR has canonical form. 
The Sasaki-Mukhanov variable C 

-C = ^' + ^Q, (3.10) 
is the gauge-invariant curvature perturbation on uniform-density hypersurfaces. From ( 13.101 1, we can express Q in terms of tp and 



(3.11) 

iC + ^)-^{C + i^) , etc. 



After plugging the above relations into ( 13.9b . and using the background equations of motion, we get a new action for the two 
variables (C, ijj): 



S2 [C, ^] = y dtd'x {ccC^ + fc CC + CAC + (^C)' + ^^C ^CA^C + ^ 

+ fip tptp + ipAip + dtp Ai/jA^V + (Ai/')^ + ^^tp ipAip + m^, ip'^ 

+fci' Ci' + /cv C'0 + CAi/' + ioci, CAV" + w^C V^AC + dc^^ACAip + d^^ ACA^V + i^Ci^ CV'} > 

(3.12) 

where the various coefficients can be found in Appendix IA 21 

It is important to note that there is no ijP' term in ( 13.12b . Moreover, the coefficient of the kinetic term for is oc (pg, and thus, 
C also becomes non-dynamical in the absence of matter field ip. In the presence of matter field there is only one dynamical 
degree of freedom in our system, which we can identify as C,. 

Eq. ( 13.121 ) is rather complicated. However, we can further simplify the expression. First we note that all the coefficients in 
( 13. 12b are evaluated in terms of the background fields, and thus are functions of time only. Then, for a general function F{t) and 
spacetime field 0, we have the following convenient relation 



j dtd^xa^F{t)4>4>:^ J dtd^x^ ("Y^'^^'*) J ^^'^^^"^ 

where "~" denotes up to total derivative terms. Similarly, we have (noticing that A = /a^) 



4 {F + 3HF 



J dtd^xa^F{t)^A(j)~ j dtd^ 



- [F + HF 



(j)A<j). 



Thus, by using the above relations and performing many integrations by parts, we find that the second-order action ( 13.121 ) can be 
recast into a rather convenient form: 

52[C>] = j dtd3xa3{c^-c' + cr4(A)c + ^ri(A)v+(r2(A)c + r3(A)c)^} , (3.13) 



where we have defined 



ri(A) = -]^{h^+ Hh^^ A + d^, A^ + di, A^ + A + - i (/^ + iHf^ 

T2{A) = (tJv< + t^C^) A + rfcV'A^ + dcV' A^ + (mc^ " ki> " SiJ/cv-) , 
r3(A) = /icvA, 

r4(A) = coc A + A^ + A3 +m^-\ (j^ + SHfA , 



(3.14) 
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for simplicity (in T3, we have used the fact /^^ = J^^, see Appendix I A 2 1 for details). Note that since A = /a^ and since 
several of the coefficients in ( |3.12| i are time-dependent, the F's are in general time-dependent. The above results should be 
understood in Fourier space, where we identify A = — /c^/ a^. 

The important point is that now ip has no time-derivatives and acts as a new constraint. The equation of motion for 1/; is 



2ri(A)i^ + r2(A)c + r3(A)c = o, 



(3.15) 



from which we can solve Tp explicitly to get 



r2(A)c + r3(A)c 
2ri(A) 



(3.16) 



After plugging ( 13.16b into (13.13b . and after some straightforward calculations, we obtain an effective quadratic action for a 
single variable 



S2[C] = / dtd 



r2 

4ri ' ^ 



r4- 



4ri 



1 d 
iJ^dt 



3 



(3.17) 



C. Subtleties of the dynamical degrees of freedom 

The important lesson to be drawn from the previous subsection is that the extra degree of freedom which appears in HL gravity 
due to the loss of space-dependent time reparametrizations as a symmetry of the theory is non-dynamical in cosmology. Hence, 
the strong coupling problem discussed recently in IitIi is not present in linear cosmological perturbation theory. Let us discuss 
this result, a result which holds for any value of A. 

The symmetry group of Hofava-Lifshitz gravity is that of "foliation-preserving" diffeomorphisms, which is smaller than the 
full diffeomorphism group of GR. Thus, one may naively expect that with smaller symmetry, we are left with less gauge artifacts 
and more physical modes. In particular, it was argued that there was an additional scalar dynamical degree of freedom in Hofava 
theory, arising from the gravity sector itself. This sounds very different from the situation in GR, where gravity has only two 
physical degrees from freedom — the two polarization states of gravitational waves. According to the standard treatment of 
cosmological perturbation theory in GR, the quantum perturbations of the scalar sector are expected to be generated by scalar 
matter fields. In this case, if there is no matter, all perturbation modes are gauge artifacts. Thus, if there is indeed one additional 
scalar dynamical degree of freedom from gravity itself in Hofava theory, our traditional picture of cosmological perturbation 
would not apply. In fact, the extra dynamical degree of freedom would lead to serious problems for the theory, as stressed in 
iIitIi . However, in this work, a detailed study of the perturbation theory shows that there is no additional dynamical scalar degree 
of freedom at all, when perturbations in the presence of matter are found by consistently expanding around the FRW metric and 
not around the flat Minkoswki one. 

Furthermore, another point is that, in Hof ava's original formulation of the theory, the lapse function N was restricted to being a 
function of time only N = N{t), where the corresponding Hamiltonian becomes non-local. In this work, we relax this restriction. 
By assuming N = N{t, x) and starting from the most general expansion of the action, we see that the apparent additional degree 
of freedom is non-dynamical. 

Is there a new dynamical degree of freedom in the absence of matter? Let us take the limit of our equations when H and the 
energy density of matter tend to zero. Since the coefficient of the kinetic term Q"^ is proportional to the background scalar-field 
value LfjQ (see Appendix IA 2b . c,^ cx i^q. Thus if there is no matter field, all perturbation modes become non-dynamical, which is 
exactly the case in GR. However, the limit we are discussing is singular unless A = 1 because if A 7^ 1 one of the coefficients in 
the action (the d^, coefficient) blows up. 

Thus, our work also shows that in the case A = 1, the case in which the action of HL gravity reduces to that of GR in the IR 
limit, there are no extra gravitational degrees of freedom in the vacuum. This result is in agreement with the conclusions in 13, 
where a scalar degree of freedom in the gravity sector was identified for A ^ 1, and where it was shown that this mode becomes 
non-dynamical when A = 1. Our conclusions also agree with the recent analysis of (0] ^. The extra scalar gravitational degree 
of freedom in the vacuum sector of the theory was also discussed in lfl6ll . That paper shows that there is a singularity for A = 1 
which is agreement with our conclusion that for this choice of A there are no new gravitational degrees of freedom. 



IV. COSMOLOGICAL PERTURBATIONS AND THE POWER SPECTRUM 

In |[l2ll (see also ifisll ) it was pointed out that a scalar field with Hofava-Lifshitz form (12.6b will obtain a scale-invariant spectrum 
of cosmological fluctuations, and it was then argued that a scale-invariant spectrum of curvature fluctuations may similarly result 



' The absence of an extra dynamical degree of freedom in the vacuum sector of HL gravity is also observed in f24l[23l . 
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independent of the equation of state of the background . This result follows from the fact that the action contains terms with six 
space derivatives. These terms dominate in the UV and yield a scale-invariant spectrum. 

In cosmology we are interested in the power spectrum of the induced curvature fluctuations. The formalism we established 
in the previous section now allows us to calculate this spectrum. We will show that for the gravitational potential ( I2.4l i of 
detailed balance form, the terms with six spatial derivatives cancel in the action for cosmological fluctuations. This happens 
independently of the detailed form of the coefficients of the higher derivative terms in the matter action. The reason for this result 
is that the terms in the gravitational action with six spatial derivatives (the terms involving the constant w) do not enter the action 
for cosmological perturbations. If we add terms to the action consistent with power counting renormalizability which break the 
detailed balance condition, terms with six spatial derivatives in the action for cosmological perturbations will survive. 

If we keep the detailed balance condition on the potential, then the leading terms in the UV contain four spatial derivatives. 
Thus, the initial power spectrum of curvature fluctuations is not scale-invariant ^. Thus, scale-invariance of the late time power 
spectrum will only arise for specific background evolutions, e.g. for inflationary expansion (as discussed later in this section) or 
for a matter bounce ifioll background. If we drop the detailed balance condition, then an initially scale-invariant curvature power 
spectrum results. 



A. Equation of motion 

The action (|3.17t for a gravitational potential satisfying the detailed balance condition has the general structure: 

S2[C] = I dtd^xa^ (-/(^ -nC 



(4.1) 



with 



7 

-n 



^ 3 

4ri 

^ 2 



^d_ /o^TsTa 
4a3 dt \ Ti 



(4.2) 



where the F's are defined in (I3.14l i. In order to write the action ( 14.1b in canonical form, we introduce the new variable 

u = a^C . 

After changing to conformal time rj (which is defined by dt — adrj) we have 



52 [C] = j di^d' 



27 



i_ 

27 
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(4.3) 



(4.4) 



where T-L = a' /a, and a prime indicates the derivative with respect to conformal time. Note that the above result should be 
understood in momentum space, that is we must make the replacement rj(A) T i{—k'^ / a?) . 
The classical equation of motion for the canonically-normalized variable u is simply 



Uk+'^^iv,k) Ufc = 0, 



with 



7 



X 

27 



X 

27 



(4.5) 



(4.6) 



where 7 and ft are defined in i4.2i and the variously introduced parameters can be found in Appendix lA2l We emphasize that 
in deriving the equation of motion ( |4.5l l-( |4~6l ), no approximation (beyond the restriction to linear perturbation theory) has been 
made and that thus the equation is exact. One can use ( I4.5b -( l4.6b as the starting point of a detailed investigation of the spectrum 
of scalar metric perturbations in Hofava-Lifshitz theory. 

In the following, we shall first investigate the UV limit of the above equations in a general background, and then study the 
evolution in an inflationary background. 



The possibility of obtaining a scale-invariant spectrum of metric perturbations in pure HL gravity was discussed in (5^. 
' As already mentioned in [121, the scale-invariance of an initial matter entropy field spectrum can induce scale-invariance of the curvature fluctuations via the 
"curvaton" mechanism 12711. 
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B. UV Limit of Perturbation Theory 

To study the perturbations in the UV limit, i.e. for k oo, we consider the terms with highest power of A in the equation of 
motion for the fluctuations. First, we note that the leading terms in the F's are 

Ti ~ A3 , r2 ~ dc^A^ , ,T3^h^^A, r4 ~ d^A^ , 7 ~ cc n (^-d^ + A3 , (4.7) 

where the ~ sign means that we are writing down only the leading term in the UV limit (but with the exact value of its coefficient). 
Inserting the coefficients (see Appendix lA2l ): 

dc = -jp- , dci, = , di, = -j^, (4.8) 



we find that 



n ~ 0. (4.9) 



In other words, the term in il, which could be naively expected to arise in Hofava-Lifshitz theory based on the terms in the 
Lagrangian, vanishes, exactly. The leading order contribution in powers of k starts with a fc^ term. As shown in lfl2ll . it is 
the k^ term which can naturally produce a scale invariant spectrum, while the k'^ term cannot. Thus, for a potential satisfiyng 
the detailed balance condition it is not true that in Hofava-Lifshitz theory a scale invariant spectrum can be produced in any 
background. 

The above cancellation of the k^ term is not an accident. One can prove that even for more general high order derivative 
terms in the scalar field Lagrangian, the above cancellation happens. To show that, we take the high order derivative part of the 
Lagrangian to be '° 

grnnA™^A'V . (4.10) 

m—l,n>7n 

The corresponding piece of the matter Lagrangian quadratic in Q is 

■ 2 

^ g™„A"QA"Q ^ ^ ^!^A™+"(^ + ^)2. (4.11) 



m— l,n>m r/i— l,n>??i 



Note that here ip is not a dynamical field, and needs to be solved for by minimizing the action. In the k ^ oo limit, there 
are no other A^ or higher spatial derivative terms, because the highest derivative in the gravity sector is A^. Thus, Eq. ( 14.111 ) 
corresponds to the dominant term in the action. Minimizing the action, the solution is 

V~-C. (4.12) 

Thus, the A^ or higher spatial derivative terms in the scalar field sector are canceled and thus do not contribute to cosmic 
perturbations. 

If one assumes that there are terms such as (pA^cp instead of A(pA^^(p in the scalar field Lagrangian, then the above proof is no 
longer valid. It is then possible to generate scale invariant perturbations for a general background in this case. However, such a 
scalar field Lagrangian breaks the shift symmetry in the gradient term, and is thus not conventionally used in the literature. 

At this stage, it is also important to recall that k^ terms in the action for cosmological perturbations will survive if one adds 
terms to the gravitational action which do not conform to the detailed balance condition. 



C. Perturbation Theory in an Inflationary Background 

Now we turn to the perturbation equation of motion i4.5i . Obviously, ( 14. 5t reduces to the standard case in IR. Thus, in 
this work, in order to investigate the possible differences of cosmological perturbation in Hofava theory from those in General 
Relativity, we focus on the UV region. For simplicity, we consider an inflationary background as usual, where the Hubble 
parameter is approximately constant. 

As discussed in the previous subsection, in the UV limit, the ~ k^ term exactly cancels out. The equation of motion takes the 
form 

< + fc^fc^ + ^ - — + APaA Uk^O. (4.13) 



See also for a study of scalar field Lifshitz actions. 
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where 

^2 _ (A - 1)^2^2 [i6i72(i _ 3^)2 ^ ^ip^2^2 ^ 2k2(A - 1)(3A - 
" " 32(3A- l)3<p2 



(4.14) 



is the analogue of the "speed of sound" and is the effective mass-square term. The expressions for both and A/^ can 
be found in Appendix lA3l We should keep in mind that ( 14.131 1 only describes the behavior of perturbation in UV limit, that is, 
the left-hand-side of (14.13b should be compensated with terms of order ©(l/fc^), 0{l/k'^), ■ ■ ■ etc., which we neglect in the 
following analysis. 

It is interesting to point out a connection with the "trans-Planckian problem" for fluctuations in inflationary cosmology. In 
1I29I1 . it was argued that the predictions of inflationary cosmology for the spectrum of cosmological perturbations are sensitive 
to hidden assumptions about physics on trans-Planckian scales. To demonstrate this point, fluctuations obeying a modified 
dispersion relation very much like the one we have obtained above was studied in ^Of]. However, in llsoll . the modified dispersion 
relation was not derived from any action principle but simply postulated. We have shown here that HL gravity yields a specific 
modified dispersion relation for scalar metric fluctuations. We will return to this point in a followup paper Isill . 

Note that in the case A = 1, the case in which the IR limit of HL gravity yields the Einstein action, S = and thus the ^ fc^ 
term in the equation of motion for the curvature fluctuations also vanishes. Thus, in this case, the case of most practical interest, 
the equation of motion for the curvature fluctuations reduces to the same form as in GR. The only difference is that the effective 
mass M takes a different from. 

In the case A 7^ 1, then in order to guarantee the stability of the perturbations in UV limit, 'E? has to be positive. If we assume 
1/3 < A < 1, from ( |4.14| i, it is easy to show this gives a "lower bound" for |i^|: 

16g^(3A-l)^+^W ,4,,. 

or, in terms of Vq: 

16^2(1 - 3A)2 + k^aV^ - 8k2(-1 + A)Vb < . (4.16) 

In the case A = 1, the stability condition on the solutions in the UV requires that the coefficient gi in the matter Lagrangian be 
negative, the sign we expect (since it is the gi term which dominates in the IR and must give the standard kinetic term for the 
scalar field Lagrangian). 

As an application of the framework we have developed in this paper, let us consider the evolution of fluctuations in an infla- 
tionary background" with a(?y) = — 1/ [Hrj), where rj is the conformal time. The equation of motion (14. 13b has been extensively 
investigated. On super-Hubble scales the solution is 

Ukiri) ^ a{-q) (4.17) 

which corresponds to constant curvature fluctuation C,. In the special case A = 1, the solution Ukiil) will be oscillating on 
sub-Hubble scales. 

More generally, the equation can be solved analytically to get'^ 

Ukiv) = ^ e-^/^V^ {^CskiiY U{a, y+l,z), (4.18) 



where C/(a, + 1, z) is a confluent hypergeometric function of the first kind with 

1 ,>2 

a=-{v + l) 



V \ 



9 a72^ (4.19) 

4 ^ 



and the overall normahzation constant is 



1 



T{a-v)T{a*) T{a)T{a*-v) 



(4.20) 



We neglect slow-roll corrections to the scale factor evolution. 

Equation of motion with the same structure of j4. 13) has also been analyzed in |33|| . in the investigation the statistical anisotropy and large-scale CMB 
anomalies. 
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where denotes a complex conjugate. Here the fc-independent normalization constant C is chosen so that Uk [rj] is normalized 
in the usual way (unit Wronskian): 

Uk{v)<{il)-ul{j^)u'^{r^)^i, (4.21) 

which is the condition for canonical quantization. 

The evolution of a perturbation mode is shown in FiglT] It can be seen that the behavior of perturbation in Hofava theory is 
very similar to that in GR. In particular, after the wavelength exits the sound horizon, the perturbation modes are frozen. On 
super-horizon scales, since 

Uk{Tl) ' r^^, , , (4.22) 



rCV(a) \iHEh] 

the dimensionless power spectrum of the curvature fluctuation ^ can easily be calculated 

Pc(fc)^|^IC(fc)P 



2^ 



Ukiv) ^ 



(4.23) 



277/ c3cc|C|r(a)r(a*) \HE 

where we have used the fact that in UV limit, 7 w C(^. The dimensionless power spectrum P({k) is shown in Fig|2](up to an 
overall factor). The figure also shows for comparison the power spectrum obtained with the standard GR mode function 

MV) = ^e'(''+^)?V^H(i)(-c,fcr?) , (4.24) 

which describes the fluctuation of a massive light scalar field in de Sitter space-time. We observe that the perturbations are 
suppressed in Hofava-Lifshitz gravity compared with those in GR. This verifies the argument that renormalizability of gravity 
generally reduces the amplitude of perturbations ll32ll . It would be also interesting to see whether this suppression of perturbations 
could stabilize the inflationary background in the UV limit and prohibit eternal inflation. 

We also would like to mention that we have not proved whether ( isa conserved quantity beyond the slow roll approximation 
here. Although in the solution, we can see ( is no more than slowly varying on super Hubble scales. In one Hubble time, the 
variation of ( is at most as large as the slow roll parameters. It is still possible for C to receive 0(1) corrections during about 60 
e-folds of inflation. Another possible contribution to C on super Hubble scales is the UV-IR transition. As the UV-IR transition 
happens at 0(1) e-folds, it may also yield a 0(1) contribution to 

Applications of the e qua tion of motion for cosmological fluctuations in HL gravity to non-inflationary backgrounds will be 
left to a followup paper IBlll . 



D. IR Limit 



Now we turn to the IR limit behavior of perturbations as an important consistency check of perturbation theory in Hofava 
gravity. In the IR Hmit (k 0), ( I4.5b -( l4.6b takes the form 



clk^ - — + M^a^ ) Ufc = , (4.25) 

a 



where the expressions for and can be found in Appendix I A 41 There are additional terms on the left-hand-side of the 
equation which are of the order 0{k'^), O(fc^), etc. which we can neglect in IR limit. Then, Eq. ( 14.251 ) has the same form as the 
corresponding equation in standard perturbation theory in GR. 

As we know, in the IR limit, Hofava theory reduces to GR with the following parameters; 



c= — -\ , WttG = — -\ , Agr = — r-^- r. (4.26) 

4 V 1 - 3A ' 8 V 1 - 3A ' 32(1 - 3A) 

Thus, setting the speed of light c in the IR to c = 1 corresponds to choosing A = ^^''^■i^l'^'' in Hofava theory. With this value of 
A, from ( |A39l l it is easy to show that 



5l = -2gi , 



(4.27) 
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qqU I ' 

0.1 1 10 100 1000 10'' 

scale factor a 

FIG. 1: Evolution of the mode function with the scale factor a. The red curve corresponds to 
the mode function in Hofava theory, the black curve shows the standard mode function in GR. 
The dashed vertical line denotes the sound horizon. The parameters were chosen to be Cj, = 1, 
k = H = 10, M = 0,3 = 0.1. 



80 
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0^ ^ ^ ' ' ' ' 

0.1 0.5 1.0 5.0 lO.O 50.0 iOO.O 

Comoving wavenumber k 

FIG. 2: Power spectrum on super-Hubble scales. The red line denotes the spectrum of perturbations 
in Hofava theory, and the black line denotes the standard GR result. The parameters were chosen 
to be Cs = 1, H = 0.1, V = 1.52, H = 10. The spectra are evaluated at the conformal time 
77 = -0.001. 

which is what we expected. Thus, the equation of motion for perturbations in Hofava-Lifshitz theory in the IR limit indeed 
reduces to the same form as in GR. However, the effective mass M and the parameter 7 still remains different from GR. One can 

check that when one further takes A = 1, we have 7 = and M is suppressed by slow roll parameters. So in the g = —1/2 
case in the IR limit, the perturbation theory completely returns to GR up to the leading order of slow roll approximation. This is 
a consistency check for our calculation. 

To summarize, we have found several cases in which the perturbation equation takes the same form as that in GR: (1) In the 
UV limit with A = 1. In this case, the equation of motion has the same form with that in GR. However, the coefficients are 
different. (2) In the IR limit with arbitrary A. In this case, the equation of motion is the same with that in GR up to leading order 
in slow roll approximation. (3) In the IR limit with A 1. In this case, the equation of motion should completely reduced to GR. 

V. CONCLUSIONS AND DISCUSSION 

In this paper we have studied the theory of linearized cosmological perturbations in Hofava-Lifshitz (HL) gravity. In this study, 
it is important to expand about a dynamical background since the presence of matter implies that the average energy density does 
not vanish. We have found that the extra degree of freedom which could be expected to arise because of the reduced symmetry 
of HL gravity is in fact not dynamical. This conclusion holds for any value of A. Taking the flat space-time limit of our analysis, 
we can also show the absence of any new dynamical scalar metric degree of freedom in the absence of matter. Our limiting 
procedure, however, only works in the special case A = 1, the case in which the HL gravitational action flows to the action of 
General Relativity (GR) in the infrared. We thus do not see any evidence of the "strong coupling problem" mentioned in UtI] . 

Starting from the most general metric including scalar cosmological perturbations, we have worked out the quadratic action for 
cosmological perturbations. It turns out that the distinguished dynamical variable for fluctuations is the usual Sasaki-Mukhanov 
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variable. In terms of this variable, the kinetic terms in the action has canonical form. The second metric variable enters the action 
without kinetic term and is hence not a dynamical degree of freedom. 

Based on the action for cosmological perturbations, we can compute the spectrum of these fluctuations. Rather surprisingly it 
turns out that the terms in the equation of motion for these fluctuations containing six spatial derivatives vanish if the potential 
term in the HL gravitational action is of detailed balance form. It is the presence of these terms which leads to the scale-invariance 
of the spectrum of HL matter lfl2l[l3ll . Thus, we conclude that the spectrum of cosmological perturbations is not scale-invariant 
in HL cosmology (with potential satisfying the detailed balance form), unlike the spectrum of spectator matter field fluctuations. 
However, the fact that spectator HL matter fields acquire a scale-invariant spectrum will likely make it possible to use the curvaton 
mechanism IztIi to induce scale-invariant fluctuations in HL cosmology independent of the expansion rate of space. 

If the gravitational action is not of detailed balance form (and this seems to be the preferred case JH [l3l if the IR limit of the 
theory is really to reproduce Einstein gravity), then the terms in the action for cosmological fluctuations will persist, making 
it possible to have a scale-invariant initial spectrum of adiabatic cosmological perturbations along the fines suggested by lfT2ll . 

We would like to warn the reader that Hofava-Lifshitz gravity faces many challenges before it can be declared as a viable 
candidate theory for quantum gravity (see lfl7l[34ll for some potential problems). We have only addressed one of these problems 
- the strong coupling problem for additional fluctuation modes, a problem which has been considered fatal for the theory - and 
shown that in fact does not arise. Other potential problems remain to be resolved. If they can be successfully resolved, then it 
becomes of great interest to explore applications of the equations of cosmological perturbations which we have derived here to 
non-inflationary backgrounds. Work on this issue is in progress lisTll . 



Acknowledgments 

We would fike to thank L. Alvarez-Gaume, G. Dvali, E. Witten and K. Zarembo for useful discussions. R.B. thanks the CERN 
Theory Division for financial support and hospitality in a wonderfully stimulating atmosphere. He also thanks C. Charmousis and 
A. Padilla for correspondence. XG and YW are supported in part by a NSFC grant No. 10535060/ A050207, a NSFC group grant 
No. 10821504, and a 973 project grant No. 2007CB815401. The research of R.B. is supported in part by a NSERC Discovery 
Grant and by funds from the Canada Research Chair program. 



APPENDIX A: VARIOUS COEFFICIENTS 
1. Coefficients in <3.9t 



8H^{- 


4 + 3A) + k2(-H 










8iJ2(- 


-1 + 3A) + k2(_14 






8ir(-l + 3A)^o 





(Al) 



(A2) 



(A3) 



K^^o{iH-3HX)^o-{-l + X)V') 
8H^i-l + 3X) + K,^i~l + X)ipl ' ^ ' 



/^^ = -Sifl , (A6) 
~ _ {-l + 3X){~n^^l + 8HV') 

J^^- SH^~l + 3X) + K^i-l + X)^y ^ ' 



hj, = 



8i?2(-l + 3A) + K2(_i + A)^; 



(A8) 
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16i?2(i - 3A)2 + 2k2(-1 + A)(-l + 3X)ipl ' 



K^A^^ (g(-l + 3A)v?o + (-1 + A)yO 
2(-l + 3A) (8i/2(-l + 3A) + + X)ipl) 



+ A)/^2 (I6ij2(^ _ 3^)2 _^ ^4^2^2 ^ (i _ 4A + 3A2) (^2) 

= 8(1 - 3A)2 {8H^~1 + 3A) + .«2(_i + x)ipl) ' ^^^^^ 



AtVo + 877^2(1 - 3X)ipoV' - 4k2(_i + A) {V'f + 32H^{1 - 3X)V" - At^ipl {3H^{-1 + 3A) + (-1 + X)V" 



8 (8i?2(-i + 3A) + K^-l + X)(pI) 



3(l2i72(i_3A) + K2(/)2) 



K2A/i2 



2H - 6HX 



^^(-1 + A)Am2^2 

2iJ(-l + 3A) (8iJ2(-l + 3A) + k2(_i + A)vj2) ' 



9wl 



(A13) 



m^^ = 3V' . (A15) 
2. Coefficients in OJD 



^ 4(-l + 3A)^2 

8iJ2(-l + 3A) + A^2(_l + A)^2' 



_ 12g(l - 3A) 3^.2 ^0^' 

^ = — ^5 — + ^^^^^ 

/C0 ^ ' (A20) 



(A21) 



hi' = .rr, . . o^^.n.U. . , o T ^ . , ^^ .... > (^22) 



(A23) 



(A24) 
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* 2 

4 = ^ , (A25) 

_ 92^1 + A)/.^ {l&H\l i\f + /.^A^A^^ + (1 - 4A + 3A^) y^g) 

m 8(1-3A)2(8H2(_i + 3A) + k2(_i + ^)^2) ' ^^^^^ 



dc^ = , (A27) 



dc^ = , (A28) 



^ 2H^nH-l + 3A)Am2 + {K^A^,^ - 8(1 - 3A)2gi) 
8(ff-3i7A)2 

= — , (A30) 

^^C^-^, (A31) 

_ (K%M^^8(l-3A)V)v^g 
'^^^ = 8{H-3H\r ' ^^^^^ 

('-3iJ«:2(^4 _^ 24F2(_i + 3A)(^o^' - 2AtVo^' + 4iJ(-l + 3A) {V'f + AH{-1 + 3X)ipl {9H^ - V"] 
= -^^ 



8iJ3(-l + 3A) 



(A33) 



aW^K (-1 + 3Aj (A34) 
-4i7(l - 3A) (sGH^il - 3A) + (V'f'j - AHk^{~1 + 3X)ipl {QH^ + V") | , 

(-kVo^' + 2i?(-l + 3A) {V'f + 2if(l - 3A)^2y//^ 



2iJ3(-l + 3A) 



3. cl and in <4J3t 



(A35) 



The exact form for the "effective speed of sound" is given by 
1 



2048i/2(l _ 3A)6g3(^4 



+ A)V^ (I6i72(i - 3A)2 + K^A 
16i7*K^ (1 - 4A + 3A2)%4 + ij2^io^_^ ^ ^ 3A)AV' 



2(^2 



+32(-l + 3A)3.93 {iQE'^K^il - 3X)'^Afi^ + (k^A/i^ _ 8(1 - 3A)2.gi) ipl {8H'^{-1 + 3A) + ^^(-i + X)ipl))] } . 

(A36) 

In the case A = 1, 

, MH^t^Kn^ + 16i/2 (k4Au2 - 32oi) ifl 
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which reduces further to = — 2gi if we set A = 0, as expected. 
The effective mass-square term AP is given by 



25oc^a^ ^ ^ ^ ^ ^ ^ 

"^d^dcip ~ '^dc^^d^d^^d^ + 2d^d^ (^Hh^d^.^ + d^ (d^jp + {^'^Hh^^ + 2 (o;^^ + 
2 (^m^rf^^ - 2d^{mQ^dQ^ + 3Hf^d,i, - 2mQd^ - iHd(^^fQ^) 



(A38) 



- 32c^ 



4. and in fOst 



5? EE 



2c^ - '^rn^) ' 



(A39) 



and 



(A40) 



with parameters are given in Appendix IA 21 
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